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Abstract 

We use the resonances of the spherical shell potential to present a thorough description of the 
Gamow (quasinormal) states within the rigged Hilbert space. It will be concluded that the natural 
setting for the Gamow states is a rigged Hilbert space whose test functions fall off at infinity faster 
than Gaussians. 
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I. INTRODUCTION 



Resonances are intrinsic properties of a quantum system, and they describe the system's 
preferred ways of decaying. The experimental fingerprints of a resonance are either a sharp 
peak in the cross section or the exponential decay of the probability to find the unstable 
particle. The sharp peaks in the cross section are characterized by the energy En at which 
they occur and by their width Fr. Decay is characterized by the energy E-^ of the particle 
and by its lifetime tr. 

The Gamow states are the wave functions of resonances, and they are eigenvectors of 
the Hamiltonian with a complex eigenvalue. The real part of the complex eigenvalue is 
associated with the energy of the resonance, and the imaginary part is associated with the 
width. The time evolution of the Gamow eigenfunctions abides by the exponential decay 
law. 

The Gamow states are able to describe both sharp peaks in the cross section and decay, 
in accordance with the phenomenological perception that resonances and unstable particles 
are two sides of the same phenomenon. As well, when the (complex) resonance energy tends 
to a (real) bound-state energy, the Gamow eigenfunction becomes a bound state, in accor- 
dance with the phenomenological perception that unstable states are only quantitatively, 
not qualitatively, different from bound states, the only difference being that unstable states 
have a non-zero width, whereas the width of stable states is zero. 

In a way, the Gamow states complete the so-called Heisenberg program, according to 
which spectral lines, widths and lifetimes are all observable quantities, and quantum me- 
chanics should be able to predict them. 

Gamow introduced the energy eigenfunction with complex eigenvalue in his paper on a- 
decay of atom, uuclei Q. and|ts propert.es and appl.cations have been considered by many 
authors, see for example |2h86| and references therein. A pedestrian introduction to these 



states can be found in 



35!]. Gamow's treatment does not fit within the Hilbert space 



though, because self-adjoint operators on a Hilbert space can only have real eigenvalues. 
Recall however that Dirac's bra-ket formalism does not fit within the Hilbert space but 
rather within the rigged Hilbert space. Similarly, the rigged Hilbert space mathematics 
asserts the legitimacy of Gamow's proposition. In the rigged Hilbert space language, the 
Gamow states are eigenvectors of the dual extension of the self-adjoint Hamiltonian. Such 
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extension can surely have complex eigenvalues 87 1. 

A rigged Hilbert space (also called a Gelfand triplet) is a triad of spaces 



$C-HC$^ (1.1) 

such that 7/ is a Hilbert space, $ is a dense subspace of and is the anti-dual space 
of The space $ has a topology that is finer than the topology inherited from "H. The 
space contains the continuous, antilinear functionals over $. Associated with the rigged 
Hilbert space (II. ip . there is always another rigged Hilbert space, 

*CHC$', (1.2) 

where is called the dual space of $ and contains the continuous, linear functionals over 
Since the space is bigger than "H^ = "H, and since is also bigger than "H' = "H, some 
physically meaningful states that find no accommodation in 1-L will find accommodation in 
and For example, the eigensolutions of the time-independent Schrodinger equation 
associated with either the scattering energies (the Lippmann-Schwinger kets \E^)) or with 
the resonant energies (the Gamow kets I^r)) find accommodation in whereas the bras 
((^E\ and (zr| find accommodation in 

The present paper is devoted to show how the rigged Hilbert space is able to accommodate 
the Gamow states. Throughout the paper, rather than working in a general setting, we will 
use the example of the spherical shell potential. 



r(x) = V{r) 



< r < a 

^0 a<r <b (1.3) 
6 < r < oo , 



and restrict ourselves to the s partial wave. However, as explained in Appendix A of Ref. |63|, 
the result is valid for any partial wave and for spherically symmetric potentials that fall off 
faster than exponentials. 

For the potential (11. 3p . expressions such as those for the Gamow eigenfunctions and the 
S matrix depend on the square root of the energy rather than on the energy itself. It is, 
therefore, easier to do calculations with the wave number k, 



k=\l^-^E, (1.4) 
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rather than with the energy E. However, we will write most results in terms of the energy, 
because they tend to be simpler than in terms of the wave number. Also, when the energy 
and the wave number become complex, we will denote them by z and q, 



2m , , 

q=\j^z, (1.5) 

and when they correspond to a resonance R, we will denote them by z-^ and /cr. 



2m 

kK=\j -^z^ . (1.6) 

We will re-write most expressions in Dirac's bra-ket notation, because of its simplicity, clarity 
and beauty. 

Wave functions in the position representation, denoted by and Gamow bras and kets, 
denoted by {z-^\ and I^r), will have sometimes a superscript + or — attached to them, and it 
is important to understand what this superscript means. Let us suppose that is a Gaussian 
wave packet in the position representation. Such Gaussian could be either an "in" state, in 
which case we denote it by y?"'", or an "out" state, in which case we denote it by y?". When 
we write the Gaussian will be expanded by the "in" Lippmann-Schwinger bras and kets, 
and its energy representation will always be the one associated with the "in" bras and kets. 
When we write the Gaussian will be expanded by the "out" Lippmann-Schwinger bras 
and kets, and its energy representation will always be the one associated with the "out" 
bras and kets. Thus, the superscripts ± are sort of "phase-space" labels, since they tell us 
which energy representation we are using, even though we may be working in the position 
representation. Physically, the superscripts ± are a reminder of whether we have imposed 
the "in" or the "out" boundary conditions on the Gaussian packet. For the Gamow bras 
and kets, the meaning of the superscripts ± is analogous. 

Of all the previous attempts to describe the Gamow states within the rigged Hilbert 



space, our approach is closest to that of Bollini et al. 23|, |24 [.There are, however, two main 



differences between the present approach and that of Refs. 



23 



2J]. First, the test functions 



we are going to use fall off at infinity faster than Gaussians, whereas the test functions 



used in Refs. 



23 



24j fall off at infinity faster than exponentials. We use test functions 



that fall off faster than Gaussians because they enable us to perform resonance expansions 



that include all the resonances of the system and that exhibit time asymmetry 



. And 



second, we obtain the relation between the Breit-Wigner amplitude and the Gamow states by 



transforming to the energy representation, whereas in Refs. 23|, |2J| such relation is obtained 
by transforming to the momentum representation. 

In Sec. ini the Gamow states of the spherical shell potential will be constructed. The 
Gamow kets associated with resonances and anti-resonances will be defined as the solutions 
of homogeneous integral equations of the Lippmann-Schwinger type. We will solve these 
integral equations in the radial, position representation. In this representation, those integral 
equations are equivalent to the time-independent Schrodinger equation subject to a purely 
outgoing boundary condition (POBC). We will also obtain the "left" Gamow eigenfunctions 
and will comment on the analogy between bound and resonance states. 

In Sec. mil we will apply the theory of distributions to construct the Gamow bras and 
kets, which in Sec. HV] will be shown to be generalized eigenvectors of the Hamiltonian 
with complex eigenvalues. Also in Sec. llVt we will construct the rigged Hilbert spaces that 
accommodate the Gamow bras and kets. 

Next, in Sec. |Vl we will obtain the energy representations of the Gamow bras and kets, 
and show that they can be written in terms of the complex delta function and the residue 
distribution. 

In Sec. IVIt we will let the energy run over the full real line in order to obtain the "energy 
representation" associated with the Breit-Wigner distribution. We will show how the com- 
plex delta function becomes the Breit-Wigner distribution in such "energy representation." 
The results of Sees. |V] and VTl will, in particular, provide a mathematical support for the 
results presented in Ref. [63 1. 

The time evolution of the Gamow bras and kets will be calculated in Sec. IVIII We will 
argue, although not fully prove, that the time evolution of a resonance ket is valid for positive 
times only, whereas the time evolution of an anti-resonance ket is valid for negative times 
only. Thus, the time evolution of resonances is given by (non-unitary) semigroups, which 
express the time asymmetry built into a decaying process. This time asymmetry seems to 



be what some authors such as Fonda et al. 



Cohen- Tannoudji et al. [90j, or Goldberger 



and Watson [91| have called the irreversibility of a decaying process. 

For the sake of completeness, in Sec. IVIIII we will construct the resonant expansions and 
see how such expansions allow us to isolate each resonance's contribution and to interpret 
the deviations from exponential decay [92 ]. 

In Sec. lIXt we will present two analogies that help to understand the physical meaning 
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of the Gamow states. The first analogy is that between the resonance expansions, the Dirac 
expansions, and the classical Fourier expansions. The second analogy is that between the 
classical, quasinormal modes and the quantum mechanical resonances. We will also explain 
the physical reason why the Gamow eigenfunctions blow up exponentially at infinity. 



II. THE GAMOW EIGENFUNCTIONS 



The Gamow eigenfunctions are customarily defined as eigensolutions of the Schrodinger 
equation subject to the POBC. Although we could start the study of the Gamow states 
with that definition, we will follow instead a treatment parallel to that of the Lippmann- 



Schwinger equation 93-951. We will define a Gamow state as the solution of an integral 

nn 

equation [11|, [1^ that has the POBC built into it. Needless to say, in the end the explicit 
solutions of that integral equation will be found by solving the Schrodinger equation subject 
to the POBC. 



A. The integral equation of the Gamow states 

The Gamow states are solutions of a homogeneous integral equation of the Lippmann- 
Schwinger type. If = -Er — irR/2 denotes the complex energy associated with a resonance 



of energy E^i and width Fr, then the corresponding Gamow state |zr) fulfills 111. Il4| 



The +iO in Eq. (12.11) means that we are working with the retarded free Green function, which 
has a purely outgoing boundary condition built into it. The retarded free Green function 
is analytically continued across the cut into the lower half plane of the second sheet of the 
Riemann surface, where the complex number z-^ is located. Therefore, as pointed out in jl^ . 
Eq. (12. ip should be written as 

This notation intends to express that we first have to calculate the retarded free Green 
function {E — Hq + iO)~^ in the physical sheet, and then continue it across the cut into the 
lower half plane of the second sheet. 
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The integral equation (12. ip has the POBC built into it. To be more precise, in the position 
representation Eq. (12. ip is equivalent to the time-independent Schrodinger equation subject 
to the condition that far away from the potential region, the solution behave as a purely 
outgoing wave. 

As is well known, to each resonance energy z^i there corresponds an anti-resonance energy 
that lies in the upper half plane of the second sheet. The integral equation satisfied by 
the anti- resonance state \z^) reads as 

1^^) = I ^n ^'^^) = i™* F H ^n ^'^) • (2.3) 
z^ — Hq — lU E^'^R E — Hq — 1\J 

In contrast to Eq. (12. ip . Eq. (12. 3p has a purely incoming boundary condition built into 
it. That is, in the position representation, Eq. (12. Sp is equivalent to the time- independent 
Schrodinger equation subject to the condition that far away from the potential region, the 
solution behave as a purely incoming wave. 



B. The Gamow states in the position representation 

In the radial, position representation, Eqs. (12. ip and (12. 3p become 

(r^R) = {r\- ^l—r^V\z^) = lim (r| ^^-^ ^1^) ' (2-4) 

(r|4) = (r| 1- = lim (r| ^ V\E) . (2.5) 

In [l^ , these integral equations are written as 

POO 

M(r;zR)=lim / G^{r, s; E)V{s)u{s; E) ds , (2.6) 

E^ZR Jo 

POO 

M(r;4)=lim / {r, s; E)V{s)u{s; E) ds , (2.7) 
Jo 

where 

u{r;zn) = {r\zn). (2.8) 

In order to obtain the explicit expressions of the Gamow eigenf unctions, instead of solv- 
ing the integral equations (12.60 and (12. 7p . we solve the equivalent Schrodinger differential 



equation 

d2 



2m dr^ 



+ V{r)]u{r;zR) = z^u{r;zR) , (2.9) 
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subject to the boundary conditions built into those integral equations, 



u{0; zn) = 0, 

u{r; zji) is continuous at r = a, 6 , 
d . 

—uir; j IS continuous at r = a, o . 
dr 

u{r; zji) ~ e'*^'^'' as r — cxd , 



(2.10) 
(2.11) 
(2.12) 
(2.13) 



where condition f l2.13p is the POBC. In Eqs. f l2.9p -f l2rT3|) . u{r; z^) = {r\z^) can denote either 
a resonance or an anti-resonance state. 



For the spherical shell potential (11. Sp . the only possible eigenvalues of Eq. (12. 9p subject 
to (I2.10p -( 12.13p are the solutions of the following transcendental equation: 



(2.14) 



where J7+ is the Jost function, see, e.g., Refs. 35|, |93|. The solutions of this equation come 
as a denumerable number of complex conjugate pairs Zn, z^. The number Zn = En — ir„/2 
is the nth resonance energy. The number z^ = En + i^nf^ is the nth anti-resonance energy. 
The corresponding resonance and anti-resonance wave numbers are given by 



_ 2m 
-k„ — \l ^^z^ , n — 1,2, 



(2.15) 



which belong, respectively, to the fourth and third quadrants of the fc-plane. For the poten- 
tial ( II. sp . the resonance poles are simple (see [3l| for an example of a potential that produces 
double poles). 

In terms of the wave number the nth Gamow eigensolution reads 



sin(fc„r) 



< r < a 



u{r; Zn) = u{r; A;„) = Nn < 



^lMgiQ„r ^ ^iMg-iQ„r a<r<b 



(2.16) 



6 < r < oo , 



where 



A'^ is a normalization factor. 



iV^ = ires [S{q)] 



q=kn 1 



(2.17) 



(2.18) 
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and Ji-Ji are coefficients whose expressions follow from the matching conditions (12. lip and 
fl2.12p . The Gamow eigensolution associated with the rath anti-resonance pole reads 



u{r- z*J = u{r; -k*J = M„ < 



j3(lfc*) sin(-/i;» < r < a 

#1?^e-'«"^' + fJ^e'Q"'^ a<r<b (2.19) 



where M„ is a normalization factor, 

M2 = ires[^(g)]^_,, =(iV2r, (2.20) 

and where 

'Qn = ^^{^n-Vo) . (2.21) 

For the sake of brevity, we will label the anti-resonance wave numbers —A;* and —Q^, the 
energies , the normalization factors M„ and the eigenfunctions u{r; z^) with a negative 
integer n as 

kn, Qn, Zn, Nn, u{r;zn) u = -1, -2, . . . . (2.22) 

This notation will enable us to write results that are true for both resonances and anti- 
resonances just once. 

Since they are eigenfunctions of a linear differential operator, the Gamow eigenfunc- 
tions (I2.16P and (I2.19P are defined up to a normalization factor. The normalization we 
have adopted was introduced by Zeldovich jsl, who used a Gaussian regulator to damp the 
exponential blowup of the Gamow eigenfunctions and obtain a meaningful normalization: 

/•oo 

lim / dre-'''''[u{r;zn)]^ = 1, (2.23) 

where n = ±1,±2.... Zeldovich's normalization has (at least) three advantages. First, 
it generalizes the normalization of bound states; second, the residue of the propagator at 
the resonance energy factors out as a product of two Gamow eigenfunctions, see Eq. f l2.29p 



below; and third, Zeldovich's normalization makes u{r; Zn) have dimensions of 1 / A/length, 
so \u{r; Znjl"^ has dimensions of a radial probability density, just like any normalized wave 
function in the position representation. 

It is worthwhile noting that the expressions for the delta-normalized Lippmann-Schwinger 
eigenfunctions are different when expressed in terms of k from when expressed in terms of 
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E QJ, |95|. However, similarly to bound states, the expressions for the normalized Gamow 
eigenf unctions are the same when expressed in terms of /cr, as when expressed in terms of 
zr, see Eqs. ^LWf and (^JU^. 



C. The "left" Gamow eigenfunctions 

After having obtained the "right" Gamow eigenfunctions, which will be associated with 
the Gamow kets, it is easy to obtain the "left" Gamow eigenfunctions, which will be asso- 
ciated with the Gamow bras. 

The "left" Gamow eigenfunctions can be obtained by complex Hermitian conjugation of 
the "right" Gamow eigenfunctions [97], or by analytic continuation of the "left" Lippmann- 



Schwinger eigenfunctions 



95 



961]. The resulting "left" Gamow eigenfunction associated with 



the resonance (or anti-resonance) energy z„ is given by 

{zn\r) = [u{r;z:)]*, n = ±l,±2,.... (2.24) 

Thus, contrary to naive expectations, the "left" Gamow eigenfunction is not just the complex 
conjugate of the "right" eigenfunction, but the complex conjugated eigenfunction evaluated 
at the complex conjugated energy. Note that this procedure to obtain the "left" from 
the "right" eigenfunctions generalizes the procedure to obtain the "left" from the "right" 
eigenfunctions of both the bound and the scattering eigenfunctions. Because the Gamow 
eigenfunctions satisfy 

[u{r-z:)]* =u{r-z^), n = ±l,±2,..., (2.25) 

the "left" and the "right" Gamow eigenfunctions are actually the same eigenfunction, 

{zn\r) = [u{r; z*J]* = uir; Zn) = {r\zn) , n = ±1, ±2, . . . . (2.26) 

In terms of the wave number, Eq. ( 12.260 reads as 

{kn\r) = [u{r;-kl)Y = u{r;kn) = {r\K), n = ±l,±2,.... (2.27) 

Note that Eq. f l2.25p is a symmetry of the Gamow eigenfunctions, and it is such symmetry 
what in the end makes the "left" eigenfunction be equal to the "right" one. Note also that 
such symmetry does in general not hold when we change the normalization of the Gamow 
eigenfunctions — yet another reason to choose Zeldovich's normalization. 
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Equation fl2.26p makes it clear why Zeldovich's normalization for the Gamow states is 
written as in fl2.23p rather than as 

POO 

lim / dre-^''''\u{r;zn)\'^ = 1. (2.28) 

Also, Eq. fl2.26p can be used to show that at a resonance (or anti-resonance) pole, the residue 
of the Green function is given by 

res [G{r, s; z)]^^^^ = —kn res [G(r, s; g)]g=fc^ = u{r] k^) u{s] k^) , n = ±1, ±2, . . . , (2.29) 

which in bra-ket notation becomes 

1 

res[(r|— -^|s)]2=^„ = (r|2;„)(z„|s) , n = ±l,±2,.... (2.30) 

Note that this factorization could have been used to define the above normalization of the 
Gamow states and to show that the "left" Gamow eigenfunction {zn\s) is the same as the 
"right" Gamow eigenfunction (s|z„). 



D. Bound states 



For the sake of simplicity in the expressions, we have chosen a potential that doesn't bind 
bound states. We would nevertheless like to briefly comment on what happens when bound 
states appear. 

The bound states satisfy the same integral equation as the resonance states, and therefore 
they automatically follow from the Schrodinger equation subject to the POBC along with 
resonances. Thus, the eigenfunction u{r; z-^i) becomes a bound state when we substitute 
the complex resonance energy by a real bound-state energy E^. In addition, Zeldovich's 
normalization for the Gamow eigenfunctions reduces to the standard normalization of bound 
states when we substitute by E^. 



III. THE GAMOW BRAS AND KETS 



The Gamow eigenfunctions u{r; z^) are obviously not square integrable, i.e., they do not 
belong to the Hilbert space L^([0, oo), dr). Thus, like the Lippmann-Schwinger eigenfunc- 



tions 



93 



95[ |. the Gamow eigenfunctions must be treated as distributions. By treating them 



as distributions, we will be able to generate the Gamow bras and kets. 
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According to the theory of distributions 



98 



eigenfunction u{r; Zn) must be defined as [OSHOS] 



the Gamow ket associated with the 



(3.1) 



LP I — {(p\zn) ■■= dr [f{r)]* u{r; Zn) , n = ±1, ±2, . . . . 

The elements (p{r) of $cxp are such that their "nice behavior" compensates the "bad be- 
havior" of u{r; Zn) so the integral f l3.ip makes sense. The space $exp will be constructed in 
Sec. IIVI In the bra-ket notation, definition fl3.ip becomes 



dr {Lp\r){r\zn) ■ 



(3.2) 



Similarly, the Gamow bras associated with the resonance (or anti-resonance) energy z„ 
are defined as 



exp 



c 



ip I — > {zn\^) := dr (p{r)u{r; Zn) , n = ±1, ±2, 



that is. 



dr {zn\r){r\ip) . 



(3.3) 



(3.4) 



From the above definitions and from Eq. f l2.26p . it follows that the actions of the Gamow 
bras and kets are related by 



(^|z„) = «|^)*, n = ±l,±2. 



(3.5) 



Since the Gamow eigenfunctions are the same when expressed in terms of the energy as 
when expressed in terms of the wave number, the Gamow bras and kets, unlike the delta- 
normalized Lippmann-Schwinger bras and kets, are the same when expressed in terms of 
the energy as when expressed in terms of the wave number: 



\kn) — l^n) ; {kn\ — (^n| ; ^ — il; i2. 



(3.6) 



IV. THE RIGGED HILBERT SPACES FOR THE GAMOW BRAS AND KETS 

Likewise any bra or ket, the Gamow bras and kets are dealt with by means of the rigged 
Hilbert space rather than just by the Hilbert space. The rigged Hilbert space we will use is 
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very similar to, although not the same as the rigged Hilbert space of Refs. 23|, |2J]. We will 
denote the rigged Hilbert space for the bras by 



$expCL2([0,oo),dr)c$^ 



exp 



(4.1) 



and the one for the kets by 



$expCL2([0,oo),dr)c$,% 



(4.2) 
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The procedure to construct the space of test functions $exp has been explained in 
95|. The most important property one has to look at is the "bad behavior" of the Gamow 
eigenfunctions. Such "bad behavior" must be compensated by the "nice behavior" of the 
elements of $exp so the integrals ( I3.ip -( P^ converge. Since the regular solution xij"^ l) of 
the Schrodinger equation is related to the Gamow eigenfunction by 

1 J^{kn 



2i Nr, 



-u{r;kn), n = ±l,±2,... 



(4.3) 



and since by, for example, Eq. (12.6) in Ref. 99| the regular solution satisfies 



|x(r;g)|<C 





q 


r 


1 + 


Q 


r 



JIm(g)|r 



qeC 



(4.4) 



the "bad behavior" of the Gamow eigenfunctions is given by 



\u{r-kr,)\<C 







\kn 


r 


\J-{kn)\ 1+1 


h 


r 



^|Im(fc„)|r 



n 



±1,±2, 



(4.5) 



Because the bound (14. 4p is sharp [99|, so is the bound (14. 5p . Thus, the Gamow eigenfunc- 
tions grow exponentially as r tends to infinity, and, in order for the integrals (I3.ip - (l3.4p to 
converge, the wave functions of $exp must fall off at infinity sufficiently rapidly. 



From Eq. (14. 5p . it is clear that the integrals in Eqs. (13. 1| 



functions that fall off at infinity faster than any exponential 



-(13. 4p converge already for 



23 



24| . Thus, exponential 



falloff is the weakest falloff that we need to require from the wave functions of $exp, see 



Refs. 



23 



24| . However, we are going to impose a stronger, Gaussian falloff because it allows 
us to perform certain resonance expansions, as will be discussed in Sec. IVIIII 

Using the estimate (14. 5p . and following the procedure of |93l-l95| to construct spaces of 
test functions, one ends up finding that $exp is given by 



^exp = {'^ eV \ \\(p\\m,ra' < OO , 171, Ul' = 0, 1 , 2, . . .} 
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(4.6) 



where T> is the maximal invariant subspace of the Hamiltonian, 

oo 

p = Pi , (4.7) 



m=0 

and II ■ ||m,m' is given by 



W\\m,m' ■= \ / dr 



OO 



— — e-'-V2(i + i7)mV(r) 
1 + mr 



m,m' = 0,1,2,... . (4.8) 



Hence, $exp is just the space of square integrable functions which belong to the maximal 
invariant subspace of H and for which the quantities (14. 8 p are finite. In particular, because 
(p{r) satisfies the estimates (14. 8p . ip{r) falls off at infinity faster than e~^^, that is, its tails 
fall off faster than Gaussians. 

Note that we have arrived at the same space of test functions as the one for the analytically 
continued Lippmann-Schwinger bras and kets [95|], since also in that case we have to tame 
real exponentials. 

Once we have constructed the space $exp, we can construct its dual ^^^p ^^d antidual 
^exp spaces as the spaces of, respectively, linear and antilinear continuous functionals over 
$exp; and therewith the rigged Hilbert spaces (14. ip and (14. 2p . The Gamow bras and kets 
are, respectively, linear and antilinear continuous functionals over $exp- As well, they are 
(generalized) eigenvectors of the Hamiltonian. 

The following proposition, whose proof follows exactly the same steps as the proof of 



Proposition 2 in [95|, encapsulates the results of this section: 



Proposition 1. The triplets of spaces ( [^^.ip and Iji4-^ rigged Hilbert spaces, and they 
satisfy all the requirements to accommodate the Gamow bras and kets. More specifically, 

(i) The II ■ ||m,m' '^'^e norms, and they define a countably normed topology, i.e., a meaning 
of sequence convergence. 

(a) The space $exp is dense in L^([0, oo), dr). 

(Hi) The space $exp is invariant under the action of the Hamiltonian, and H is $exp- 
continuous. 



(iv) The kets \zn) are continuous, antilinear functionals er $exp) i-^-, \zn) ^ 



exp 
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(v) The kets \zn) are generalized "right" eigenvectors of H with eigenvalue 

H\Zn) = Zn\Zn) , 71 = ±1, ±2, . . . ; (4.9) 

that is, 

{(p\H\Zn) , ip G $exp • (4-10) 

(vi) The bras {zn\ are continuous, linear Junctionals over ^^xp, i-^-, {zn\ € ^exp- 

(vii) The bras {zn\ are generalized "left" eigenvectors of H with eigenvalue Zn'- 

{Zn\H = Zn{Zn\ , n = ±1, ±2, . . . ; (4.11) 

that is. 

If), y^e^cxp. (4.12) 

Proposition 1 makes it clear, in particular, that there is a 1:1 correspondence between 
Gamow bras and kets. 

Note that in terms of the wave number, the eigenequations (14. 9 p and (14.111) become 

H\K) = ^kl\K), n = ±l,±2,..., (4.13) 

{kr,\H = ^kl{kr,\, n = ±l,±2,.... (4.14) 
Note also that the bra eigenequation (14. lip is not given by 

{Zn\H = zl{Zn\ , (4.15) 

as one may naively obtain by Hermitian conjugation of the ket eigenequation (14. 9p . The 
reason lies in that one has to use complex Hermitian conjugation to obtain the "left" from 
the "right" Gamow eigenfunction, see Eq. (I2.26p . 

The normalization condition satisfied by the Gamow states is the following: 

{Zn\Zn>) = 5n,n' , n, n' = ±1, ±2, . . . . (4.16) 

When n = n', Eq. (I4.16P follows from Zeldovich's regularization (I2.23p . When n ^ n' , 
Eq. (I4.16P can be proved in the same way as one proves the orthogonality of bound states: 

{Zn\H\Zn') = Zn{Zn\Zn') = Zn'{Zn\Zn') , (4.17) 
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where we have made use of the fact that {zn\ and \zni) are eigenvectors of H with eigenvalue 
Zn and Zn', respectively. From the second equahty in fl4.17p . we obtain 

{Zn - Zn'){Zn\Zn') = 0, (4.18) 

which yields the desired result, since z„ 7^ Zn' when n ^ n! . It should be noted however that, 
similar to the normalization of scattering states, the normalization condition f l4.16p has only 
formal meaning and does not imply the use of a Hilbert-space scalar product. For example, 
the "scalar product" built on Eqs. (14.161) and (12.231) would not satisfy (/, /) > 0. 



V. THE ENERGY REPRESENTATIONS OF RIGGED HILBERT SPACES AND 
OF THE GAMOW BRAS AND KETS 

We turn now to obtain and characterize the energy representations of the rigged Hilbert 
spaces (14. ip and (14. 2 p and of the Gamow bras and kets. It is here where we will need to 
introduce the labels ± in the notation for the wave functions and for the Gamow bras and 
kets. 



A. The energy representations of the rigged Hilbert space 

The "in" and the "out" energy representations of $exp are readily obtained by means of 



the unitary operators IJ± of 



9^: 



t^i^exp = $±exp , (5.1) 

which in turn yield the energy representations of the rigged Hilbert spaces (14. ip and (14. 2p : 

$±exp C L2([0, 00), d^) C $'±,,p , (5.2) 

8±exp C L2([0, 00), dE) C $Lxp • (5.3) 

The elements of $±exp will be denoted by (p^i^z) = U±(p{z). 

In 95|, we characterized the analytic and growth properties of the wave functions in the 
wave number representations, (p^{q), which are related to the wave functions in the energy 
representations as 
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Thus, the results of 95| also characterize the analytic and growth properties of (p^{z), and 



we will refer to 



95| whenever we need to make use of any such properties. 



As mentioned above, from now on we will add a label to the action of the Gamow states. 



n 



±1,±2 



(5.5) 



When we use the label +, it will mean that the energy representation is obtained through 
the operator f/+, and when we use the label — , it will mean that the energy representation 
is obtained through the operator f/_. 



B. The energy representations of the Gamow bras and kets 

In order to obtain the energy representations of the Gamow bras and kets, we first need 
to define the linear complex delta functional at z: 

{5.\ : *exp ^ C ^^^^^ 

^ I — > (41^) := ^{z) , 

where $exp may be either $+oxp or $_exp, and (p may be either cp^ or Lp~ . Thus, the linear 
complex delta functional at z associates with each test function, the value of the test function 
at z. One can write (15. 6 p as an integral operator as 

(S,\ip) = / dESiE - zME) = (piz) . (5.7) 
Jo 

In this way, one can interpret the complex delta function 6{E—z) as the analytic continuation 
of the Dirac delta function 6{E — E'). 

The antilinear complex delta functional \Sz) at the complex number z can be defined in 
a similar way: 

1^.) : ^exp ^ C 

^ ^ := mz*)Y . 

We also need to define the linear and antilinear residue functionals at z: 



(5.8) 



(res^l : $exp i — > 



(p I — )■ (fes^l^) := res[^(2;)] 

|fesz) : ^exp I — > C 

^ I — )■ (^|fes4 := ies[(p{z*)Y 



(5.9) 



(5.10) 
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where res[^(z)] stands for the residue of (p at z. Both the complex delta functionals and 
the residue functionals at z are well defined when the test functions can be analytically 
continued into 2, as is our case jS]. 

We will need also the following normalization factor: 

Ml = iTes[S{z)U,„ = i|!-2fc„res[^(g)],=fc„ = ^2KN'^ , (5.11) 

where was used in Sec. [TTl to normalize the Gamow eigenfunctions. 
If we denote the energy representations of the Gamow bras and kets as 

{^Zn\ = (^zj?7±, (5.12) 

|2^) = t/±kn), (5.13) 
then the following proposition, whose proof can be found in Appendix \^ holds: 

Proposition 2. For a resonance (or anti-resonance) of energy Zn, the "minus" (or "out") 
energy representation of the Gamow bras and kets is given by 

(res.J, n = ±l,±2,... , (5.14) 



■/^ n 

|2;) = iy27K|5,„), n = ±l,±2,.... (5.15) 
Their "plus" energy representation is given by 



-Zn\=lV2^Mn{5,^, n = ±l,±2,..., (5.16) 
72^ 



l^) = -^|f^s,J, n = ±l,±2,.... (5.17) 



Proposition 2 shows, in particular, that the "plus" energy representation of a Gamow bra 
or ket is different from its "minus" energy representation, thereby showing that the labels 
± matter. 

The complex delta functional and the residue functional can be written in more familiar 
terms as follows. By using the resolutions of the identity 



/■oo 

/ = / dE\E^){^E\, (5.1^ 
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we can formally write the actions of ( integral operators and obtain 



poo 

/ dE{^Zr,\E^){^E\cp^) 

Jo 

POO 

/ dE{^z^\E^)(p^iE). (5.19) 
Jo 







Comparison of fl5.19p with fl5.14p and fl5.16p shows that ( Zn\E ) is proportional to the 
residue distribution, 

{-Zn\E-) = -^ies[-l„, E>0, n = ±l,±2,..., (5.20) 

and that is proportional to the complex delta function, 

{^Zn\E^) =iV2^Afn5iE - Zr,) , E>0, n = ±l,±2,.... (5.21) 

Similarly, by using fl5.18p we can formally write the actions of as integral operators: 

{^^\zi) = {y.^\zi) 

^dE{cp^\E^){^E\zi) 

POO 

= / dE[^^{E)r{^E\zi). (5.22) 

By comparing (15.220 with (I5.15P and (15.171) . we deduce that {~E\z~) is proportional to the 
complex delta function 

{-E\z-) = iV2^MnS{E - Zn) , E>0, n = ±l,±2,... , (5.23) 

and that ("^'-E'l^;^) is proportional to the residue distribution, 

(+E|z+) = -^res[-],„, E>0, n = ±l,±2,.... (5.24) 

It is important to realize that with a given test function, the complex delta function 
and the residue distribution at Zn associate, respectively, the value and the residue of the 
analytic continuation of the test function at z„. This is why when those distributions act on 
[(p^{E)]* as in Eq. fl5.22p . the final result is respectively ["^"(-2*)]* and res [^+(2;*)]*, rather 
than [i^~(2;„)]* and res [^+(z„)]*, since the analytic continuation of [(p^{E)]* is [i^(z*)]* 
rather than [^(2;)]*. 
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VI. THE (-00,00)- "ENERGY" REPRESENTATION 



The spectrum of our Hamiltonian is [0, 00). Hence, in Eqs. fl5.2ip and fl5.23p the energy 
E runs over the positive real hne. In this section, we are going to let E run over the full real 
line. In doing so, we can see what would happen if the spectrum of the Hamiltonian wasn't 
bounded from below. 

It is important to keep in mind that in this section, we will need to treat resonances 
and anti-resonances separately. Also, strictly speaking, whenever we say that E runs over 
over the full real line (—00,00), it will actually mean that in the case of resonances (anti- 
resonances), E runs infinitesimally below (above) the real axis of the second sheet of the 
Riemann surface. 

In order to construct the (—00, 00)- "energy" representation, we will construct the trans- 
form A that lets the energy vary over the full real line. The transform A is modeled after 



the "6' transform" of 



19|, and it allows us to connect the physical spectrum, which in our 



case coincides with [0, 00), with the support of the Breit-Wigner amplitude, which coincides 
with (—00, 00). Basically, A takes a test function Lp^{E), E > 0, of $±cxp into its analytic 
continuation over the full real line, (p^{E), E G (—00, 00). In order to distinguish when the 
energy runs over the physical spectrum from when it runs over the full real line, we will 
denote ip^{E), E G (—00, 00), by '^^{E) and thus will write 

A^^ = ^^, (6.1) 

and 

A$±exp = ^±cxp = {V^{E) I E G (-00, 00)} . (6.2) 

The following diagram shows how A links the energy representation with the (—00, 00)- 
"energy" representation: 

$±exp C L^([0, 00), di?) C $±cxp energy representation 
i A ; A (6.3) 

$±exp C L2(M,d„E) C $±cxp (-00, 00)- "energy" repr. 
where L^(M, d^i?) is the following space: 

L2(M,d„E) = {/ I /GL2([0,oo),dE), lim / dE |e-'^"/(E)|2 < 00} . (6.4) 

a— >0 



00 
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In Eq. (16 ■4p . the integral is assumed to be calculated in the second sheet, infinitesimally below 
(or above, in the case of anti-resonances) the real axis. The convergence factor e~'^" (which 
becomes e^"^"" in the case of anti-resonances) is needed because the analytic continuation of 
(p^{E) into the negative energies blows up exponentially [95]. Actually, if it wasn't needed, 
the spectrum would be the full real line. Nevertheless, the space L^(M, do-E) is not crucial 
to our discussion. 

It is important to understand that although we have denoted the functions ^ and 
ip^ = A(p^ by a different symbol, they are indeed the same function. More precisely, 
they are different "pieces" of the same function. In particular, the value of their analytic 
continuation at a complex number z is the same, 

^^{z)=^^{z). (6.5) 

Obviously, the analytic continuation of their complex conjugates enjoys an analogous prop- 
erty, 

[^^{z*)r = [^^{z*)r. (6.6) 

We use different symbols for different "pieces" of the same function because the proof of the 
connection between the Breit-Wigner amplitude and the complex delta function becomes 
more transparent. For resonances, such connection is given by 

A|2;) = |^^-), n = l,2,..., (6.7) 

JZ/ Zfi 

where the ket \^ — ~) is associated with the Breit-Wigner amplitude as follows: 



1 -\ 



I E — Zn I ^■'^P 

r ^ (^-|^-) :=lim.^ondi?e-^"(-^^) [^-(E)]^ 



(6.8) 



We will call this ket the Breit-Wigner ket. The integral in Eq. (16. 8 p is supposed to be 
calculated in the lower half plane of the second sheet, infinitesimally below the real axis. By 
the properties of ^p~{z) in the lower half plane of the second sheet 95|], the Breit-Wigner ket 
is a well defined antilinear functional. The proof of (16. 7p is provided in Appendix |Al 

The combination of (16. 7p with the results of Sec. |V] shows that the Gamow eigenfunc- 
tion u{r;Zn), the complex delta function (multiplied by a normalization factor) and the 
Breit-Wigner amplitude (multiplied by a normalization factor) are the same distribution in 
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different representations: 

u{r; Zn) ^ iV2^^J{E - z„) , E e [0, oo) ^ --^ , E G (-cx), cx)) 

(6.9) 

posit, repr. energy repr. (— oo, oo)- "energy" repr. 

Physically, these links mean that the Gamow states yield a decay amplitude given by the 
complex delta function, and that such decay amplitude can be approximated by the Breit- 
Wigner amplitude when we can ignore the lower bound of the energy, i.e., when the resonance 
is so far from the threshold that we can safely assume that the energy runs over the full real 
line. However, because there is actually a lower bound for the energy, the decay amplitude is 
never exactly given by the Breit-Wigner amplitude. Mathematically, the reason lies in that 
A is not unitary, which makes the energy representation be not equivalent to the (—00, 00)- 
"energy" representation. 

One can also relate the "plus" Gamow bra with a Breit-Wigner bra: 

,, 1 



E - z„ 

where the Breit-Wigner bra is defined as 



n = l,2,..., (6.10) 



(6.11) 



The proof of (I6.10p is almost identical to the proof of (16. 7p . 

For the anti-resonance energies, we obtain similar results to (16. 7p and (I6.10p . The Gamow 
ket of an anti-resonance is related to a Breit-Wigner ket as 

A|$;) = |^^-), n = -l,-2,..., (6.12) 

where now the ket | ^^^ ~) is associated with the Breit-Wigner amplitude as follows: 

" (6 13) 

(ri^-) :=lim„^oroodEe'^"(^^) [^-(E)]* . 

Similarly, the Gamow bra is associated with a Breit-Wigner bra as 

(+z„|A=(+-^|, n = -l,-2,..., (6.14) 

E - Zn 
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where the Breit-Wigner bra is now defined as 

(6.15) 



$+exp ^ 



The proofs of fl6.12p and (16.1 4p are very similar to those of (16. 7p and (I6.10p . In Eqs. (I6.13P 
and (I6.15p . the integration is supposed to be done infinitesimally above the real axis of the 
second sheet, contrary to Eqs. (16. 8p and (16. lip , where the integration is supposed to be done 
infinitesimally below the real axis of the second sheet. Also, in Eqs. (I6.13P and (I6.15P the 
regulator is e'^", a > 0, whereas in Eqs. (16. 8p and (16. lip the regulator is e~'^", a > 0. The 
reason why anti-resonances need the opposite sign in their regulator will become apparent 
in Sec. IVml 

Note that unlike \z^) and {'^Zn\, the "plus" Gamow ket \z^) and the "minus" Gamow 
bra (" not related to a Breit-Wigner amplitude in an obvious way. 

The relation between the various representations we have constructed can be conveniently 
summarized in diagrams. For resonances we have 

H;if'{r) $exp C L\[0,oo),dr) C {r\Zn) = u{r; Zn) 

if/, i f/_ i 

H;ip-{E) $_exp C L2([0,oo),dE) c $!!exp {-E\z-)=i^/2^K6{E-Zrr) (6.16) 

lA I A 

and 

H;ip+{r) $exp C L2([0,oo),dr) C {-^ Zn\r) = u{r; z^) 

IU+ iu+ iu+ 

H;ip+{E) $+exp C L2([0,oo),dE) C S^exp {-^Zr.\E+) = iV2^Mj{E-Zn) (6.17) 
i A I A 
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where H denotes the operator multiphcation hy E, E > 0, and H denotes the operator 
multiphcation byi?, — oo<-E'<oo. The top, middle and bottom rows of these diagrams 
contain, respectively, the position, the energy and the (— oo, oo)- "energy" representations. 
For anti-resonances, the diagrams are analogous. 



VII. THE TIME EVOLUTION OF THE GAMOW STATES 



We are now going to obtain the time evolution of the Gamow states by extending the 
time evolution operator q-^^^/^ into the spaces ^^^p and $exp- Since such extension was 



constructed in [95| to obtain the time evolution of the analytically continued Lippmann- 
Schwinger bras and kets, and since the Gamow states can be obtained from the analytically 
continued Lippmann-Schwinger bras and kets, the time evolution of the Gamow states will 



easily follow from the results of 95 1. 



Let us calculate first the time evolution of the Gamow bra {'^Zn\ for a resonance energy: 



Zr 



\Q^Ht/h^^ , t < only 



that is. 



= e^^"*/^'^(z„) , t < only 

= e^^"*/^(+z„|<^+) , t < only , G $+exp ; (7.1) 

{+Zn\e"'^''^ = e'^-'/^i^Znl , only for t < , n = 1, 2, . . . . (7.2) 



The reason why the time evolution for a Gamow bra associated with a resonant energy Zn 
is defined only for t < is that when t > and Zn = E.^ — ir„/2, the factor e'^"*/'* blows up 
exponentially, and therefore e'-^^/^^Lp^ violates the bound (18.61) below. Hence, e^^*^'^ $'^ i s not 



in $+exp when t > 0, and therefore the dual extension of e^^*/^ is not well defined lOOj. We 
should also note that, strictly speaking, Eq. (17.11) does not prove that the time evolution of 
{^Zn\ is well defined for t < in the sense of the theory of distributions. In order to prove 
so, one needs that the space of test functions ^exp be invariant under the action of e^^*/'' for 
t < 0. Since it is not known whether $exp is invariant under the action of q^^^/^^ it remains 
an open problem to show that Eq. (17. 2p holds in the sense of the theory of distributions. 
Let us now calculate the time evolution of the Gamow bra {^Zn\ for an anti-resonance 
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energy: 



that is, 



= (+2-Je^^*/V^) , t>Oonly 
= e''-'l^^{zn) , t > only 

= e''-"\''Zn\^-^) , t>0 only , G $+e.p ; (7.3) 

(+;z„|e-'^*/^ = e'""*/^(+z„| , only for t > , n = -1, -2, . . . . (7.4) 



The reason why the time evolution for a Gamow bra associated with an anti-resonant energy 
Zn = En + ir„/2 is defined only for t > is that when t < 0, the factor e^^"*''^ blows up 
exponentially, and therefore e'^*/''^"'' violates the bound (18. 6p . Hence, e^^^/^( p'^ is not in 



+exp when t < 0, and therefore the dual extension of is not well defined |lOO| . As in 
the case of Eq. (17. ip . Eq. (17. 3p does not prove that the time evolution of {^Zn\ is well defined 
for t > in the sense of the theory of distributions when Zn is an anti-resonance energy. 
The time evolution of the Gamow ket \z~) associated with a resonant energy Zn is given 

by 

(^-|e-^*/^1.-) = (^1e-^*/^1z;) 

= (e'<*/V"«))* , i > only 
= e-^^"*/''(^-«))* , t>Oonly 

= e-^^"*/^(^-|z-), t>Oonly, G $+e.p ; (7.5) 

that is, 

^~''"''^\Zn) = e-''"'/^\z-) , only for t > , n = 1, 2, . . . . (7.6) 

Similarly to Eqs. (17. 2p and (I7.4p . Eq. (17. 6p is clearly not defined for t < 0, although it 
remains to be proved that it holds for t > in a distributional way. 

When we consider an anti-resonance, it can be easily shown that Eq. (17. 6p becomes 

e-'^'/^\z-) = e-''"'/^\z-) , only for t < , n = -1, -2, . . . . (7.7) 

Similarly to Eqs. (Q, (ED and (ESD, Eq. ((TZD is clearly not defined for t > 0, although it 
remains to be proved that it holds for t < in a distributional way. 
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In summary, the time evolution of the Gamow states is given by non-unitary semigroups 
and therefore is time asymmetric, expressing the irreversibihty of a decaying process. Such 
semigroups are simply (retarded or advanced) propagators that incorporate causal boundary 



conditions through the analytical properties of the test functions 



95| . However, as explained 



above, the rigorous proof of f l7.2p . f l7.4p . (17. 6p and fl7.7p is still lacking, because it is not known 



whether $exp is invariant under e '^*/'*. 



VIII. RESONANCE EXPANSIONS 

The Lippmann-Schwinger bras and kets are basis vectors that were used to expand nor- 



malizable, smooth wave functions in 94j |: 



oo 

±\/±77l,.±\ 



{r\ip^)= dE {r\E^){^E\ip^) . (8.1) 







The Gamow states are also basis vectors. The expansion generated by the Gamow states is 
called the resonance expansion. 

A given quantity (wave function, amplitude, etc.) can be expanded by resonance states 
in many different ways, depending on how many resonances we include in the expansion, 
see e.g. review \4^. When we include only a few resonances close to the real axis, as in 
Berggren's and Berggren-like resonance expansions, the wave functions ip{r) must fall off at 



infinity faster than exponentials 
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24| . However, when we include all the resonances, we 
will see that the wave functions must fall off faster than Gaussians. 

For the sake of simplicity, we will focus on the resonance expansion of the transition 
amplitude from an "in" state ip~^ into an "out" state (f~: 

POO 

((^-,yp+)= / dE{^-\E-)S{E){+E\^+), (8.2) 
Jo 

where S{E) is the S matrix. For the spherical shell potential, and also for any spherically 
symmetric potential that falls off faster than exponentials, the S-matrix and the Lippmann- 
Schwinger eigenfunctions can be analytically continued to the whole complex plane (see 



Appendix A of Ref. [63|], and references therein). Thus, by using the contour of Fig. [T], we 
obtain 

-)=J2^^-K){^Zn\v+)+ / dE{^-\E-)S{E){+E\^+), (8.3) 
1 Jo 



n=l 
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where we have tacitly assumed that {!f~\E~)S{E){'^E\(p~^) tends to zero in the infinite arc 
of the lower half plane of the second sheet. The integral in Eq. (18. 3p is done infinitesimally 
below the negative real semiaxis of the second sheet. By omitting (p~ in (18. 3p . we obtain 
the resonance expansion of the "in" wave functions, 

^^ = J2 + / dE \E-)S{E){-^E\^+) . (8.4) 

n=l 

The resonance expansion for the "out" wave function can be obtained in a similar way. In 
Eqs. (18. 3 p and (18. 4p . the infinite sum exhibits explicitly the contribution from the resonances, 
while the integral is the non-resonant background. 

In obtaining Eqs. (18. 3p and (18. 4p . we have tacitly assumed that {(p~\E~)S{E){'^E\(p^) 
tends to zero in the infinite arc of the lower half plane of the second sheet. However, as 



shown in 



95|, (v9 \E )S(E)(~^E\lp~^) diverges exponentially there, since for any (3 > there 



is a constant C such that 



95| 



\{ip-\z-)\<C\q\-^/h^r , (8.5) 

I (+2|</'+) I < C\q\-^/\^^ , (8.6) 

where q is the corresponding complex wave number in the fourth quadrant of the fc-plane. 
Therefore, Eqs. (18.31) and (18.41) need to be established properly. In order to do so, one 
has to control the exponential blowups (18. 5p and (18. 6 p by calculating the time evolution of 
Eqs. ([83D and 

(yp-,e-'^*/V) = y"e-'^"*^^((^-|2-)(+^„|v9+) + / dEe-''''/''{ip-\E-)SiE){+E\ip+), 
1 Jo 

t > only, (8.7) 



n=l 



Ve-^^"*/^|z-)(+z„|(^+) + / dEe-'^'/'\E^)SiE){^E\^^) , 

n=l •>'0 

t > only. (8.8) 



These equations are valid because the following limits hold in the infinite arc of the lower 
half plane of the second sheet for any a > 0: 

lim e-'^'iy^-lz-) = lim e-^""(+z|<^+) = 0, (8.9) 
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which in turn follow from Eqs. (18. 5p and (18. 6 P (see however 100|). Equations (18. 3 P and (18. 4p 



should then be understood as the limit of Eqs. (18. 7p and (18. 8p when t — )■ 0+. 



As shown in 
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24l | , the Gamow bras and kets are already well defined when the tails of 
the test functions fall off like exponentials rather than like Gaussians. The reason why we 
chose a Gaussian falloff has finally become clear. For test functions with exponential falloff, 
the above resonance expansions make no sense, since there is no way w e can regularize 



the blowup of such test functions in the infinite arc of the second sheet |10lj . However, 
imposing a Gaussian falloff on the elements of $exp enables us to regularize their blowup in 
the complex energy plane by using the time evolution phase e"^^*/'' as a regulator. Also, it 
is clear that Gaussian falloff is the slowest falloff that can be regularized in this way. 

As is well known, resonance expansions allow us to understand the deviations from ex- 
ponential decay. If a particular resonance, say resonance 1, is dominant, then Eq. (18. 8p can 
be written as 

Q-iHt/h^+ ^ e-i^^i/'-^\z-){+z^\y^+) + background(l) , (8.10) 

where the term "background(l)" carries the contributions not associated with resonance 1, 
including those from other resonances. Because "background(l)" never vanishes, there are 
always deviations from exponential decay. The exponential law holds only when the wave 
function is well tuned around the Gamow state ), in which case "background(l)" can be 
neglected and only the resonance (Gamow state) contribution to the probability needs to be 
taken into account. 



IX. PHYSICAL MEANING OF THE GAMOW STATES 

We are now going to explore the physical meaning of the Gamow states. We will do 
so by way of two analogies. The first analogy is that between classical Fourier expansions, 
quantum completeness relations and resonance expansions. The second analogy is that 
between the Gamow states and the quasinormal modes of classical systems. As always when 
one draws analogies between classical and quantum mechanics, one should keep in mind 
that in classical mechanics the solutions of the wave equations are actual waves, whereas in 
quantum mechanics the solutions of the Schrodinger equation are probability amplitudes. 
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A. Plane waves, the Lippmann-Schwinger bras and kets, and the Gamow states 



Plane waves e'*^^ represent monochromatic light pulses of well-defined wave number k. 
Experimentally, one cannot prepare monochromatic plane waves: all that one can prepare 
are wave packets (p{k) that have some wave-number spread. The corresponding wave packet 
in the position representation, f{x), can be expanded in terms of the plane waves as 

^{x) = ^= [ dke'''''(p{k) , (9.1) 



/2tt 

which in Dirac's notation is written as 

= / dk{x\k){k\yD) . (9.2) 



When (p{k) is highly peaked around a particular wave number k^, the wave packet is well 
approximated by a monochromatic plane wave, ip{x) ~ q^^^ox^ 

The Lippmann-Schwinger bras and kets are a quantum version of the classical plane 
waves. The monoenergetic eigenfunctions {r\E^) represent a particle with a sharply defined 
energy E (and with additional "in" or "out" boundary conditions). In analogy to the Fourier 
expansion of wave packets in terms of classical plane waves, Eq. fl9.ip . the eigenfunctions 
{r\E^) expand wave functions (f^ as in Eq. fIS.ip . When the wave packet (p^{E) is highly 
peaked around a particular energy Eq, then the approximation '^^{r) ~ (r \E^) holds. 

The physical meaning of the Gamow states is similar. Likewise the monoenergetic scatter- 
ing states, the Gamow states cannot be prepared experimentally: All that can be prepared is 
a wave packet ip^ . In complete analogy to the expansions (19.1 p and (18. ip . the Gamow states 
and an additional set of "background" states expand a wave function see Eq. (18. 4p . 
When the wave function is finely tuned around one resonance, say resonance 1, then in 



general the approximation ^^{r) ~ {t\zi) holds for all practical purposes 102| . It is in this 
sense that a lone Gamow state is the wave function of a quantum decaying particle. 

When the approximation '^^{r) ~ {'r\zi) holds, the Gamow state can be used to char- 
ity in a time- dependent description of resonant scattering 



acterize the transport of probabi 
and decay. For instance, in Ref. 50|] Garcia- Gal deron et al. present the example of a delta- 
shell potential where one resonance dominates the decay of the system. In order to show so, 
the authors of |50| calculate the survival probability using a square integrable function ip^ . 



They also use a resonant expansion to approximate by one single resonant state (r|zi). 
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As shown in Fig. 3 of Ref. 50|, the exponential decay of the survival probability calculated 
by way of the Gamow state {r\zi) is indistinguishable from the one calculated by way of the 
"exact" square integrable wave function ip^ . 



B. Quasinormal modes vs. resonance states 

In classical mechanics, confined linear oscillating systems -e.g., finite strings, membranes 
or cavities filled with electromagnetic radiation- have preferred states of motion. Such states 
of motion are called normal modes. Each normal mode is associated with a characteristic 
real frequency. Unless it is perturbed, a system in a normal mode will keep vibrating 
the same way perpetually. When friction or dissipation enters into play and therefore the 
system dissipates energy, the system has preferred ways of doing so, which are called the 
quasinormal modes. Unconfined linear oscillating systems also have quasinormal modes, 
and they are obtained by i mpo sing Sommerfeld's radiation condition, which is the classical 



counterpart of the POBC |103l |. Each quasinormal mode is associated with a characteristic 
complex frequency, whose imaginary part is associated with the exponential damping of the 
oscillation. 

In quantum mechanics, normal modes correspond to bound states, and quasinormal 
modes correspond to resonance states. Much like quasinormal modes describe the system's 
preferred ways of dissipating energy, the Gamow states describe the system's preferred ways 
of decaying. The imaginary part of the complex energy of the Gamow state is associated 
with exponential decay, in analogy to the imaginary part of the complex, classical frequency 
being associated with exponential dissipation. 



C. Physical meaning of the exponential blowup of the Gamow states 

The Gamow states blow up exponentially at infinity, and it is important to understand the 
physical origin of such exponential blowup. Let us consider first the Lippmann-Schwinger 
eigenfunction {r\E^). These time-independent, non-normalizable eigenfunctions are inter- 
preted as an incoming plane wave that impinges on a target and an outgoing wave multiplied 
by the S matrix. However, the actual expression of {r\E^) does not lead to such interpre- 
tation. Only when one views the Lippmann-Schwinger eigenfunction in a time-dependent 
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fashion, one can arrive at such interpretation. Thus, even though they are time- independent, 
the Lippmann-Schwinger eigenfunctions encode what happens in a scattering experiment at 
all times. 

Similarly, the Gamow states are time- independent, non-normalizable eigenfunctions that 
describe the decay of a quantum system at all times. Since after a long time (formally, when 
t — )■ oo) the resonance will surely have decayed and gone to infinity, the Gamow state needs 
to provide a time-independent probability amplitude of finding the particle at infinity that 
is much greater than the probability of finding the particle anywhere else in space, hence 
the exponential blowup at infinity. 



X. CONCLUSIONS 



We have used the spherical shell potential to present a systematic procedure to construct 
the rigged Hilbert space of the Gamow states. A Gamow state has been defined as the 
solution of the homogeneous integral equation introduced in HQ . Such integral equation 
is of the Lippmann-Schwinger type, and is equivalent to the Schrodinger equation subject 
to the POBC. 

By applying the theory of distributions, we have constructed the Gamow bras and kets 
and shown that they are, respectively, linear and antilinear functionals over the space of 
test functions $exp; where the elements of $exp are smooth functions that fall off faster than 
Gaussians. We have shown that the Gamow bras and kets are, respectively, "left" and "right" 
eigenvectors of the Hamiltonian, and that their associated eigenvalues coincide with the 
resonance energies. We have argued, although not rigorously proved, that the exponential 
time evolution is given by a non-unitary semigroup. Such semigroup time evolution exhibits 
the time asymmetry of a decaying process. 

We have also constructed the energy representations of the Gamow states. We have 
shown that such energy representations are given by either the complex delta function or 
by the residue distribution. These results complement the properties of the Gamow states 



in the momentum representation obtained in 



14 
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Because in the position representation the wave functions in $exp fall off faster than Gaus- 
sians, we have been able to construct resonance expansions that include all the resonances. 
Such resonance expansions also exhibit the time asymmetry of the decaying process. 
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Finally, we have clarified some of the physical properties of the Gamow states by drawing 
analogies with classical Fourier expansions and quasinormal modes. We have also clarified 
the origin of the exponential blowup of a Gamow state at infinity. 
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Appendix A: Proofs 

Here we list the proofs of some results we stated in the paper. In the proofs, whenever 
an operator A is acting on the bras, we will use the notation A', and whenever it is acting 
on the kets, we will use the notation A^ : 



(±^„|AV^) := (^;^n|AV^) , V(/.± G $exp , (Al) 

:= (At^±|z±), e $exp. (A2) 

Thus, A' denotes the dual extension of A acting to the left on the elements of ^exp; whereas 
A^ denotes the antidual extension of A acting to the right on the elements of $exp- This no- 
tation stresses that A is acting outside the Hilbert space and specifies toward what direction 
the operator is acting, thereby making the proofs more transparent. 

Proof of Proposition 2. 

The proofs of Eqs. fl5.14l) - fl5.17p all follow the same pattern. We start by proving (15.151) . 
The Gamow eigenfunction u{r; Zn) is proportional to the analytic continuation of the 



Lippmann-Schwinger eigenfunction X [r; E) 55| , 



u{r; Zn) = i^l-nMnX {r; Zn) ■ (A3) 

From this equation and from the analytic properties of the elements ip~ G $-exp obtained 
in 



95|, it follows that 
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which proves flS.lSp . The proof of fl5.16l) is analogous. 



In order to prove fl5.17p . we need the following relation |55|: 



u{r; Zn) 



Mr,. 



res [x+(r;2;)]^=^„ . 



Then, 



dr [ip^{r)]*u{r; z., 



V2^ 

Mn JO 



by (O 

dr [(p"*"(r)]*res [x'*'(r; z)]z^ 
res [^{z*J]* 

(^|fes^„), V^+G$+exp, 



which proves fl5.17p . The proof of (15.140 is analogous. 



Proof of Eq. (I67fl) . 

Let (p~ G $-exp- It was proved in 



dr[ip-ir)]*uir;zn) hy 

PCX) 

W2^Un dr[^-{r)]*x-{r-Zn) by (jA3 

— exp ) 



by dMl) 



(A4) 



(A5) 



(A6) 



95] that for any /3 > 0, the following estimate is valid 



in the lower half plane of the second sheet: 



. |Im(g)|- 



(A7) 



where q is the corresponding complex wave number in the fourth quadrant of the fc-plane. 
This estimate implies that in the infinite arc of the lower half plane of the second sheet, the 
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following limit holds for any a > (see however [lOO]): 

lim e-'°^[^-(z*)]* = 0. (A8) 

\z\—^oo 

Then, by Cauchy's formula, 

e-'-^-r = / dEe-^'^''- [lf-{E)Y . (A9) 

Multiplying this equation by i-\/27rA/'„ yields 

e-^-"iv^K[^-«)]* = I" dEe-^"^ (~^) E^: ■ ^^^^^ 

From Eqs. (ESD, and (lATOj) it follows that 



iy2^Ar„[^-«)]* = ((^-|^^-) . (All) 
We now define the action of on \'z~) by 

(^-|AX|S;):=(A-V-|2;). (A12) 

Since 

(A-V-|^) = (^12;) = iV2^Afn[^{z:)Y = iv^Ar„[^-«)]* , (A13) 

we have that 

(^-|AX|2;) = (^-|-^-) , G $_exp, (A14) 



E - Zn 



which proves (16. 7p . 
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FIG. 1: The contour to obtain resonance expansions in the A;-plane (left) and in the i?-plane (right). 
It is assumed that the contour encloses all the resonances in the lower half plane of the second 
sheet, and that the radius of the arc is sent to infinity. The filled (hollow) dots represent the 
resonance (anti-resonance) poles. 
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